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1. Introduction 

Effective Lagrangians have become a widely used tool in understanding physics in- 
volving a mass gap in the spectrum. They can be used in theories in a weakly coupled 
regime but with unknown underlying physics (as is the case in the Higgs sector of the 
standard model) or in theories with a strongly coupled regime where the usual per- 
turbation formalism breaks down. Our interest will be focused on Quantum Chromo 
Dynamics (QCD). It is a well established theory at high energy where our theoretical 
knowledge and the experimental outcome agree with rather good accuracy. At low 
energy the situation is less satisfactory because the theory becomes strongly coupled 
and non-perturbative, standard perturbative methods can not be applied. One of 
the immediate differences is given by the degrees of freedom at low and high energy. 
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At the former they are characterised by hadrons while in the latter they are better 
described in terms of the fundamental interacting quarks and gluons. 

A suitable method to tackle phenomenology at low energy in the mesonic sector, 
besides direct numerical computation as done in lattice QCD, is to use the fact that 
QCD possesses an almost exact symmetry. One can then rely on these symmetries 
and their breaking pattern using an effective Lagrangian method. We will use the 
chiral symmetry present in the QCD Lagrangian in the limit of massless quarks. The 
use of this symmetry and the effective Lagrangian method is now known as Chiral 
Perturbation Theory (ChPT). It was introduced by Weinberg [1] and systematized 
by Gasser and Leutwyler [2] for the case of the light up and down quarks as well 
as for the case where the strange quark is treated as light in addition [3]. They 
performed a basic set of next-to-leading order (NLO) calculations allowing a first 
determination of all low-energy constants at NLO, the L[, invoking the Zweig rule to 
set L\ = Lg PS 0. It was hoped this could be tested in K e4 easily, but it turned out, 
when the explicit calculations were performed, that these only had very suppressed 
contributions to the form-factors [4, 5]. 

One question is how to order the various terms in the chiral expansions. The 
proposal used by most people is to count energies and momenta as a small parameter 
of order p and the quark masses as order p 2 . Alternative countings, taking the quark 
masses also as order p are possible, see [6] and references therein. Combining the two- 
flavour two-loop calculations of nn scattering [7, 8] and the pion scalar form-factor 
[9] with the Roy equation analysis [10] it could be shown that in the two flavour case 
the correct counting was the standard one [11, 12] using the recent determination of 
the pion scattering length from [13, 14]. 

The up and down quark masses are much smaller than the strange quark mass. 
The question thus remains whether three flavour ChPT itself converges sufficiently 
well to be of practical use and whether alternative countings of contributions involv- 
ing the strange quark mass need to be used. This possibility is discussed in the recent 
work [15, 16]. The argument is that disconnected loop contributions from strange 
quarks, via kaons and etas, can be large, making a convergent three flavour ChPT 
difficult to achieve in the usual sense [16, 17]. Answering this question is part of the 
larger problem of whether the strange quark can be considered a light quark. This 
was part of the motivation behind the recent work in three flavour ChPT at NNLO 
on 7T7T scattering [18] and the various scalar form-factors [19]. Earlier calculations 
of the pseudoscalar masses and decay constants, see Ref. [20] and references therein, 
showed the possibility of this behaviour. The various vector form-factors calculated 
did not seem to have problems with convergence [21, 22, 23, 24]. 

Work on ttK scattering began soon after Weinberg's calculation of 7T7T scattering 
[25]. The earliest reference known to us is [26]. During the 1970s there was a 
dedicated series of experiments culminating in the review by Lang [27]. These were 
used extensively together with dispersion relations and crossing symmetry in [28] 
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and [29]. 

The development of ChPT led after some time to the calculation of the 7iK scat- 
tering amplitude to NLO [30, 31]. There were also earlier attempts at unitarization 
of current algebra for this process. An example can be found in the discussion in 
[32, 33] and references therein. Other recent related works are the various attempts 
at putting in resonances in this process starting with [34]. Approaches involving 
resummations can be found in [35] and [36]. An alternative approach is to consider 
the kaon as heavy and treat only the pion as a Goldstone boson. This is known 
as heavy kaon ChPT. The applications to nK scattering can be found in [37] and 
[38]. Unfortunately this approach has many free parameters and does not allow to 
connect nK scattering to other processes. It does, however, have the possibility of 
being applicable even if standard ChPT does not work. 

On the dispersive side, the analyses of [28] were slowly updated to get at a 
determination of the p 4 LECs. The first work was [39] and the full analysis has 
recently become available [40]. In the mean time, the isospin breaking corrections to 
irK scattering at NLO have been evaluated in [41, 42, 43, 44]. Since we work in the 
isospin limit and the dispersive calculation of [40] has been performed in the same 
approximation we do not discuss these works further. 

In this paper we calculate the standard ChPT expression for nK scattering to 
next-to-next-to-leading-order. A large number of calculations to this order exist in 
three flavour ChPT and we have thus been able to determine most p 4 LECs with 
this precision after making some assumptions on the values of the p 6 constants. In 
earlier work it has been found that the Zweig rule suppressed constants L\ and L r 6 
could be sufficiently different from zero that the scenario of large corrections due 
to disconnected strange quark loops was not ruled out. Pushing this calculation to 
NNLO allows then to perform this comparison at the same footing as all the other p 4 
LECs. Earlier attempts at using had led to rather large errors for these constants, 
e.g. L\ = (—0.2 ± 0.9) 10~ 3 [45]. The work on nn scattering [18] and scalar form- 
factors [19] at NNLO order gave an indication that the region with L\ m 0.45 10~ 3 
was preferred. This fitted with the NLO nK work done earlier [40]. As discussed 
below, contrary to our expectations, the results from nK scattering at NNLO are 
more indicative of a smaller value for L\. 

This paper is organized as follows. In Sect. 2 we give a very short overview 
of ChPT and the references for NNLO techniques. In Sect. 3 we discuss a few 
general properties of the nK scattering amplitude. Sect. 4 gives an overview of 
our main result, the calculation of the nK scattering amplitude to NNLO in three 
flavour ChPT. We also present here some plots showing the importance of the various 
contributions. The inputs we use to do the numerical analysis are described in Sect. 7. 
The main numerical analysis is presented in Sect. 8 and we give our conclusions in 
Sect. 9. 
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2. Chiral Perturbation Theory 

ChPT is the effective field theory for QCD at low energies introduced in [1, 2, 3]. 
Introductory references are in Ref. [46]. The usual expansion is in quark masses and 
meson momenta generically labeled p and assumes m q ~ p 2 . The Lagrangian for the 
strong and semi leptonic mesonic sector to NNLO can be written as 

C = C 2 + C 4 + C 6 , (2.1) 

where the subscript refers to the chiral order. The lowest order Lagrangian is 

C 2 = ^f{uX + X + ). (2.2) 

The mesonic fields enter in a non-linear fashion via u = exp (iM/(F v / 2)), with M 
parametrising the pseudoscalar fields. The quantity also contains the external 
vector and axial-vector (a M ) currents 

u^ = i{u ] d^u- d^uv) - ivtr^u + iul^) , Z M (r>) = - (+)a M . (2.3) 

The scalar (s) and pseudo scalar (p) currents are contained in 

X± = u] X^ ± ux^u , x = 2B (s + ip) . (2.4) 

The p 4 or NLO Lagrangian, £4, was introduced in Ref. [3] and is of the general form 

12 

£4 = L ^ = ■ ■ ■ + (X+) + L ^X + ? + ■■■■ (2.5) 

i=i 

We have explicitly shown two terms with chiral symmetry breaking which in addition 
present a double flavour trace structure, which indicates that these two terms are 
suppressed by the Zweig rule. 

The schematic form of the NNLO Lagrangian in the three flavour case is 

£ 6 = C i°i ( 2 - 6 ) 

i=l,94 

and we refer to [47] for their explicit expressions. 

The ultra-violet divergences produced by loop diagrams of order p 4 and p & cancel 
in the process of renormalization with the divergences extracted from the low-energy 
constants Lj's and Cj's. We use dimensional regularization and the modified minimal 
subtraction (MS) version usually used in ChPT. An extensive description of the 
regularization and renormalization procedure can be found in Refs. [8] and [48]. The 
divergences are known in general [48, 2, 3, 49] and their cancellation is a check on 
our calculation. 
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3. The ttK amplitude: general properties 

The iiK scattering amplitude in isospin channel / can be written as 

A((n( Pl )K(p 2 ) - n(p 3 )K(p 4 )) = T\s^u) . (3.1) 

s, t, u are the usual Mandelstam variables 

s = (pi+P2) 2 , t=(Pi-P3) 2 and u=( Pl -p 4 ) 2 . (3.2) 

There are two possible isospin combinations 1 = 1/2 and / = 3/2. These two are 
related via s <-> u crossing which yields 

T*(s,t,u) = - l -Ti(s,t,u) + ^T'Hu,t,s) . (3.3) 

We also define the crossing symmetric amplitudes T + and as, 

T ± (s,t,u) = ±(±T^s,t,u)+Tl(u,t,sj) . (3.4) 

These amplitudes can be calculated most easily from the purely I = 3/2 process 
TI+K+ -> n + K + . 

In order to describe scattering kinematics it is convenient to introduce the vari- 
able 

2 s ( (m K + m 7T ) 2 \ f (m K -m n ) 2 \ 

The kinematical variables t, u can be expressed in terms of s and cos 6 as 

t = -2<£ K (1 -cos 9) , u = -s-t + 2m 2 K + 2ml. (3.6) 
The various amplitudes are expanded in partial waves via 

oo 

T 7 (s, t, u) = 16tt J2( 2i + l)^(cos 6)t\(s) . (3.7) 

Near threshold these can be expanded in a Taylor series 

= \V~sql e K (ai + b\ q l K + 0(q$ K )) , (3.8) 

defining the threshold parameters a\ and b\. 

Below the inelastic threshold the partial waves satisfy 

Im^( S ) = ^|^( S )| 2 . (3.9) 
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In ChPT the inelasticity only starts at order p 8 . In this regime the partial waves can 
be written in terms of the phase-shifts as 



The amplitudes are often expanded around the point t — 0, s — u via 

oo 

T + (s,t,u)= J2 c P j » 2j > 

i,j=0 
oo 

T-(s,t,u) = ^r,/V^ ; . (3.11) 

i,j=0 

where v = (s — u)/ (4m^) and the cfj are referred to as subthreshold expansion 
parameters. They are normally quoted in units of the relevant power of m^+. We 
will always list these parameters in increasing order corresponding to powers of t and 
s — u. 

The p 2 result for the amplitude has only a few nonzero items. As a consequence 
the imaginary parts for all other partial waves starts only at order p 8 . This allows 
the amplitude to be written in the form 

Tl(s,t,u) = Gi(s) + G 2 (t) + G 3 (u) + (s-u)G 4 (t) + (s-t)G 5 (u). (3.12) 

The functions Gi(s) have a polynomial ambiguity due to s + t + u = 2m 2 K + 2m 2 K . The 
functions Gi(s) have various singularities. G±, G3 and G5 contain singularities from 
the irK and r\K intermediate states and G2 and G4 from the possible nonstrange 
two meson intermediates. The precise relation with the various singularities can be 
found in [39]. 



4. ChPT results 

4.1 Results at order p 2 

The lowest order result is very simple and corresponds to 

c oo = °> 4) = and c oo = • (4- 1 ) 

7T TV 

All higher terms vanish. This was initially performed using current algebra methods 



4.2 Results at order p 4 

The next-to-leading calculation was first performed in [30, 31]. We present it here in 
a slightly different form, but the final expression given in that reference agrees with 
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ours up to terms of order p 6 . Note that as mentioned also in [39] there are some 
misprints in the formulas in [30]. 

We present the analytical expressions for the G; t functions defined in (3.12) and 
where the polynomial part was isolated in a function G e (s,t,u). We also use A = 
m 2 K — m 2 and £ = m\ + m 2 . 



F$ G 1 (s) = B(ml, m 2 K , s) (-1/2 s £ + 1/4 s 2 + 1/4E 2 ) , (4.2) 
Ft G 2 (t) = 5 (m 2 , m 2 , t) (1/16 1 A - 1/16 1 £ + 1/4 1 2 ) 
+B(m 2 K ,m 2 K ,t) (3/16 1 2 ) 

+B(m 2 ri ,m 2 ,t) (1/24 m\ A - 1/24 m 2 £ - 1/16 i A 

+ 1/16 1 E + 1/72 A E - 1 /144 A 2 - 1 /144 E 2 ) , (4.3) 
i^G 3 (u) = 5(m 2 , m^, u) (-3/8 «E+ 11/32 u 2 + 1/8AE 

-7/32A 2 + 1/8E 2 ) 
+B(m 2 K , m 2 , u) (- 1/16 m 2 u - 41 /288 m 2 A + 1 /32 m 2 E 



+ 1/96 + 5/96 uA- 3/32 -u E + 3/32 
+31/192 A E + 493/3456 A 2 + 3/128 E 2 ) 
+B 1 (ml,m 2 K ,u) (-1/4AE - 3/8A 2 ) 
+fii(m^, m 2 , u) (-1/8 m 2 A - 3/16 A E - 13/48 A 2 ) 
+B 21 (m 2 n ,m 2 K ,u) (3/8 A 2 ) 

+fi 21 (m 2 Y ,m 2 M )(3/8A 2 ), (4.4) 



"77 

i^G 4 (*) = fi(m 2 ,m 2 ,t) (-1/24* - 1/12A + 1/12E) 

+B(m|-, m^, t) (-1/48 1 + 1/24A + 1/24E) , (4.5) 
F$G b {u) = B(ml,m 2 K ,u) (-l/32u - 1/32A + 1/16E) 

+B(m 2 K , m 2 , u) (1/32 m 2 - 1/32 « + 3/64 A + 3/64 E) 
+fi 1 (m 2 ,m 2 oM )(l/16A) 

+B 1 (m 2 K , m 2 , u) (1/16 m 2 - 1/32 A - 1/32 E) , (4.6) 
^G 6 ( S ,t )M ) = F 2 (-l/2 S +l/2£) 

+ (16tt 2 )- 1 (3/32 m 2 s- 1/32 m 2 (t - u) - 1/16 m 2 £ 

-1/64 s A + 31 /192 s E + 1 /16 s 2 - 7/192 (i - u) A 
+ 1/64 (* - u) E - 1/48 (t - u) 2 + 1/96 AS- 13/96 E 2 ) 

+L\ (-4 s (t - u) + 2 s 2 + 2(t - w) 2 - 8A 2 ) 

+L r 2 (2 s (t - u) - 8 s E + 5 s 2 + (t - w) 2 + 4E 2 ) 

+L^ ( s 2 + {t - u) 2 - 2A 2 ) 

+L\ (—4 s £ + 4(t - u)E + 8A 2 ) 

+L£ (2sA-2sE - 2AE + 2A 2 ) 
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+L r 6 (-8A 2 + 8S 2 ; 
+LI (-4A 2 + 4£ 2 ^ 



+A(m 2 ) - 21/32 s - 5/96 (t - u) - 7/32 A + 13/24 S) 
+A(m 2 K ) (1/96 m 2 - 3/16 s + 1/48 (t - u) + 31/576 A + 5/64 £) 
+A(m 2 ) (9/160 m 2 + 3/32 s + 1/32 (*-«) + 7/64 A - 49/320 S) . 

(4.7) 

The finite part of the one- loop integrals are denoted by A, Bj as defined in [20]. 
4.3 Results at order p 6 

The full result at order p 6 is rather cumbersome 1 . In this section we quote only the 
dependence on the order p 6 constants. This contribution can be written exactly in the 
form of the subthreshold expansion (3.11) with as nonzero combinations, normalized 
to F£ m 2 J+ 2j and F* m^ 2j+1 for c+- and c~- respectively, 

c+ = 16 m 2 K+ mU {C[ + AC r 2 + C r 5 + C r 6 + 2 C 7 - C[ 2 - 4 C[ 3 - 2 C[ 4 - 3 C[ 5 

- 4 C[ 6 + 3 C[ 9 + 6 C r 20 + 1 2 C r 21 - C r 26 - 4 C r 28 + 3 + 6 C r 32 ) 



+16 m 2 + (C[ + 4 C 2 r + 2 C 6 r + 2 C 7 r + C r 8 - C[ 2 - 4 C[ 3 - 4 C[ 5 - 4 (7 



16 



—2 C[ 7 + 3 C[ 9 + 8 + 24 C 21 — C 26 — 4 C 2g + 3 C 31 + 8 C 32 ) , 
c+ = 8 m 2 K+ m 2 + (-4 C[ - 16 C 2 r + 2 C7J - C 5 r - 3 C 6 r - 4 C r 7 - C r 8 + 4 C[ 2 + 10 C r 



13 



+u 14 + o o 15 -t- z o 17 — z o 22 — z o 25 — o 26 — z o 29 — 4 o 30 j 
+8 (-CT - 4 C 2 r - 2 Q r - 2 C 7 r - C 8 r - 2 C[ 2 - 4 C[ 3 
-C[ 4 + 4C[ 6 + C 2 r 6 + 4C 2 g 



+8 m*+ (-C[ - 4 C 2 r - CI -CI -2 C r 7 - 2 C[ 2 

-2 C[ 3 + C[ 4 + 1/2 C[ 5 + 4 C[ 6 - C[ 7 + CI, + 4 C^) , 
c 00 = 16 m x + (C[ 5 + 2 C[ 7 ) 

+ 16m 3 x+ m3 + (A CI + 2C[ 4 + 2C[ 5 - 4C£ 2 + 4C 2 r 5 + 2C 2 r 6 -4^) , 
4o = "4+ (12 C[ + 48 C r 2 -8C r 4 + CI + 10 C 6 r + 8 C 7 r + 4 C 8 r + C[ 
+4 C r u - 2 C[ 2 - 4 C[ 3 + 2 C r 22 - 4 + 4 ) 
+m 2 + (12 C[ + 48 C 2 r - 8 C\ + 4 C 5 r + 5 C r e +8C r 7 + C r 8 + C[ 
+ 2 - 2 C[ 2 - 1 C[ 3 + 2 C 2 r 2 - 4 C 2 r 3 + 4 C r 25 ) , 
c- = 8 m K+ m 2 + (-4 C£ - C 6 r - + C[ + 2 - 2 C[ 2 - 6 C[ 3 + 2 C£ 2 - 2 CQ 
+8 (-4 CI -CI- 2 CI + C[ + 4 - 2 C[ 2 - 12 C[ 3 + 2 C 2 r 2 

-2C3b), 

c+ = 16 m 2 K+ m 2 + (C 6 r + C r 8 + C[ + 2 - 2 C[ 2 - 2 C[ 3 + 2 C 2 r 2 + 4 Cy 



x It can be obtained from the website [50] or from the authors upon request. 
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+ 16m 4 K+ ml + (CI + 2C r 6 + C[ + 4C r u - 2C{ 2 - 4C[ 3 + 2C r 22 + 4C r 23 ) , 
4o = 1/2 (-7 C[ - 32 C r 2 + 2CI + 10 C 4 r ) , 
c 20 = 6m K+ (-C[ + 2C r 3 + 2CI) , 
c+ = 8 m 2 K+ (3 C[ + 6 CI - 2 C\) , 

coi = 32 rr4 + (-C[ + 2 C 3 r + 2 C 4 r ) . (4.8) 
Notice that the combinations (— C[ + 2 Cg + 2 C 4 ) shows up in both c^ and c^. 



5. Resonance estimate of the contribution from the p b con- 
stants 

Up to now we relied only on chiral symmetries to calculate the amplitude function 
at low energy. In order to give an estimate of the p 6 LECs, we assume our process 
to be saturated by the exchange of vector and scalar meson resonances. The general 
formalism of resonance saturation (RS) in ChPT was described in [51], [52]. The 
places where comparisons with experiment are available are in general in reasonable 
agreement with the estimates obtained via RS. 

For both types of exchange, we only consider the polynomial contributions to 
7^7r-scattering starting at 0(p 6 ), thus directly corresponding to the C['s LEC's con- 
tribution. 

The vector resonances are included through the matrix of fields [53] with 
Lagrangian 

U = -\(V»„Vn + l -ml(V^) - |^(VKX]> + /x(^K,X-]) (5.1) 
while for the scalar meson nonet, the matrix of fields S, we consider 

C s = ^(VSV^S - M 2 S S 2 ) + c d (Su^) + c m (S X+ ) (5.2) 

After integration of the resonance fields, the Lagrangians relevant to the present case 
read 

Cv = ( fa V M K, «"] - f x V2 [«", x _]) ^ (5.3) 

£s = ^4 <(<*V"(«X) + c m V^ + ) 2 > (5.4) 

S 

where we use [8] 

f x = -0.025, g v = 0.09, c m = 42 MeV, c d = 32 MeV, (5.5) 
and the masses are 

m v = m p = 0.77 GeV, m A = m ai = 1.23 GeV, m s = 0.98 GeV. (5.6) 
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The numerical results from both contributions to the subthreshold expansion pa- 
rameters are listed in Table 1. The contribution to the full amplitude at order p 6 
corresponds precisely to the expansion (3.11) including only these subthreshold con- 
stants. The explicit expressions for the nonzero constants are 

c+ = -16 f x g v /(M 2 v V2) {m\ m\ + m% m 2 ) + 4 c 2 JM 4 s [m\ m\ + m 4 K m 2 n ) 

-16 f 2 /M 2 (m\ ml + m A K m 2 ) - 2 g v /M v (m\ m 4 . + m 4 K m 2 ) , 
4> = § fx 9v/(M v V2) (2 m\ ml + m A K + ml) - c d c m /M A s (12 m\ ml) 

+c 2 d /Ml(4mlml) + 2c 2 jMl(3m 2 K ml-m%-ml) 

+8 fjM v (m 2 K ml+m i K + m\) + g 2 v /M 2 (3 m\ ml + m A K + m\) , 
c - = c d c m /Ml (16 m\ ml) + c 2 JM a s (8 m\ ml) 

+f 2 /M 2 (96 m\ ml) - g 2 v /M 2 (4 m\ ml) , 
4o = -3 f x gv / {M 2 V2) (m 2 K + m 2 n ) +5/2 c d c m /Ml(m 2 K + m 2 n ) 

-7/4 c 2 d /M 4 s (m 2 K + ml) - 7/8 g 2 v /M 2 {m\ + ml) , 
c7 = 24 f x g v /(M 2 V2) {m K m\ + m 3 K ) - 4,c d c m /M 4 s {m K ml + m 3 K ) 

-2 c 2 d /M A s (m K ml + m 3 K ) + 3 g 2 v /M 2 (m K ml + m 3 K ) , 
c+ = 16 / x ^/(M^ y/2) (m 2 K m 2 w + m A K ) + 8 q c m /M| (m^ ml + m^) 

+4 c^/Mj {m\ ml + m 4 ^) + 2 gy/M 2 (m 2 K m 2 w + m A K ) , 
4 = c 2 d /Ml (7/8) + g v /M v (3/16) , 
c 2 "o = c^/M| (3/2 m K ) + g v /M v (3/4m K ) , 
c+ = -6cyMl(m 2 K )+gy/My(m 2 K ), 

c m = c 2 d /Ml (8m 3 K ) + g 2 v /M 2 (4m 3 K ) . (5.7) 

Even if many of our results only get a small contribution from the above RS 
arguments, these estimates are in general a major source of uncertainty in the 0(p 6 ) 
terms. The estimates from resonance exchange for the masses and decay constants 
and the related sigma terms are the most uncertain because of the simple treatment 
of the scalar sector. These are discussed in more detail in [20] and [54] . Here we have 
set many effects, e.g. the d m term of [20], equal to zero, the naive size estimate of [20] 
led to anomalously large NNLO corrections. The estimates of the amplitudes 
can be found in [53] after the work of [55]. The effect of varying these was studied 
in [53] and found to be reasonable. 

The above procedure is obviously subtraction point dependent and is normally 
only performed to leading order in the expansion in 1/N C , with iV c the number of 
colours. Many other approaches exist, some recent relevant papers addressing this 
issue are [54, 56, 57] and references therein. A systematic study of this issue is clearly 
important, for our present first study the estimates are sufficient. 
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Table 1: The resonance contributions to the subthreshold parameters. The units are 
rn^+ 2 ^ and m 2 ^ 2 ^ 1 for cjj and c~- respectively. We have also shown the chiral order 
at which they first have tree level contributions as well as the contributions with the 
L\ = CI = at n = 0.77 GeV. 



6. A first numerical look 

In this section we present a first look at the numerical results for the two loop 
amplitudes. We choose as input the pion decay constant, the charged pion mass, the 
charged kaon mass and the physical eta mass. 

F n = 92.4 MeV , m n = m n + = 139.56995 MeV , 
m K = m K+ = 493.677 MeV , m v = 547.3 MeV . (6.1) 

The subtraction scale \i = 770 MeV is used throughout the paper unless otherwise 
mentioned explicitly. 

We present the results for the subthreshold parameters with all L\ and C\ set 
equal to zero at the scale of the rho mass. Notice that the very different sizes of the 
various quantities are to a large extent given by their normalization in powers of 
and m n . 

The order p 4 results differ somewhat from those quoted in [30] and [39]. The 
reason for this is that some variation in the choice of precisely what is called p 4 and 
p 6 is possible. The choice we have made is different from those in the mentioned 
papers. Especially Cqq suffers from this numerically. E.g. taking the eta mass given 
by the Gell-Mann-Okubo formula instead, the numerical p 4 result for it changes to 
0.16. We conclude that we are also in numerical agreement with those papers. 

The higher order corrections look very large when looked upon as the contri- 
butions to the various terms in the expansions. This is partly due to cancellations 
making some quantities very sensitive to higher order corrections. 

We also present some plots of the S channel partial waves both in the I = 1/2 
and 1 = 3/2 channel. These are shown in respectively Fig. 1 and Fig. 2. Two points 
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loops only L[ = C[ = 




1 1 1 1 1 

0.25 0.3 0.35 0.4 0.45 0.5 

s [GeV 2 ] 

Figure 1: The corrections at order p 2 , p 4 and p e to the S partial wave in the 1=1/2 
channel with all LECs set to zero at [i = 0.77 GeV. 

loops only L[ = C[ = 




channel with all LECs set to zero at fi = 0.77 GeV. 

of interest are s = 0.263 GeV 2 for the subthreshold expansion and the threshold 
at s = 0.401 GeV 2 . From the sizes of the contributions at orders p 2 , p 4 and p 6 as 
shown it is obvious that there seems to be a better convergence near threshold than 
at the subthreshold point. Contrary to the tttt scattering case, the lowest order result 
has already nonlinearities. The amplitude at this level is perfectly linear in s, t and 
u but in taking the partial waves, t and u depend nonlinearly on s via q^ K . The 
partial waves have been extracted from the amplitudes using a five point Gaussian 
integration over cos#. 
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In order to show the conver- 
gence also around threshold we 
present as well in Table 2 the 
contributions at order p 2 , p 4 and 
p 6 the scattering lengths and ranges 
and the value of the T + ampli- 
tude at the Cheng-Dashen point 

s — u — m 2 K , t = 2m 2 , 

(6.2) 

as well. These will also be stud- 
ied in more detail later on when 
we add the contributions from 
the LECs and compare to ex- 
perimental and dispersion rela- 
tion results. 
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0. 
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T+ 
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1 


.141 
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010 





,831 


0.011 



Table 2: The contributions at order p 2 , p 4 and p 6 for 
the scattering lengths and ranges and the amplitude 
7. Input parameters at thg C heng-Dashen point with the LECs set equal 

to zero as well as the resonance estimate to order p 6 . 

For our SU (3) ChPT results we 
use as inputs the masses and de- 
cay constants given in Sect. 6 and a subtraction constant /i = 770 MeV. We work in 
the limit of exact isospin. 

In addition we use the full refit of the L\ {% = 1,2,3,5,7,8) to order p 6 us- 
ing a range of values for L\ and L\ as input. These values of L\,i = 1, . . . , 8, are 
used to evaluate the matrix element. This is the same procedure used in [19] and 
[18] to study the variation of some observables as a function of the vector (L r A) L r & ). 
The experimental inputs used for the fitting procedure are: i) the values of the 
form-factors as measured by the E865 experiment [13, 14], f s (0), g p (0), f' s (0), g' p (0), 
ii) the pseudoscalar decay constants F w ±, F K ± and Hi) the masses of the pseudoscalar 
mesons, m n ±,mK±,'m r) . The performed fits correspond to fit 10 in [45] but with dif- 
ferent input values for the vector (LJ, Lg). These represent the only free parameters, 
in the analysis. The quark-mass ratio m s /m used is fixed to be 24. The variation 
inside the range (20 — 30) was studied in [45]. 

The estimates of the p % contributions to nK scattering we use are those given 
above. These lead to the central value contributions to the various threshold pa- 
rameters given in Table 2. The uncertainty on these is quite considerable. Other 
resonance estimates of ttK scattering can be found in [34, 35]. 

The iiK scattering amplitude also obeys relations from crossing and unitarity. A 
new recent analysis using these methods is [40] . Once the choice of the high energy 
input is done, no more freedom is allowed. The constraints at the matching point 
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are stronger than in the case of the Roy analysis of ttty scattering. We also quote for 
comparison the results from the older analysis of [28]. This is what we use as our 
main "experimental" input for 7rK scattering. 



8. Numerical Analysis 

8.1 7i K only 

Let us first look at the sub- 
threshold expansions and com- 
pare the dispersive calculations 
with our results. The results from 
the two analyses [28] and [40] 
can be found in Table 3 together 
with our calculation for the L\ 
corresponding to fit 10 of [45]. 

The results for c\ Q and Cq , 
which are the only two that ob- 
tain a lowest order contribution, 
are shown in Figs. 3(a) and 3(b). 
Cq has a large lowest order con- 





Fit 10 


[40] 


[28] 


c 00 


0.278 


2.01 ± 1.10 


-0.52 ±2.03 


c 10 


0.898 


0.87 ±0.08 


0.55 ±0.07 


c oo 


8.99 


8.92 ±0.38 


7.31 ±0.90 


c 20 


0.003 


0.024 ±0.006 




c 10 


0.088 


0.31 ±0.01 


0.21 ±0.04 


c 01 


3.8 


2.07 ±0.10 


2.06 ±0.22 


c 30 


0.0025 


0.0034 ±0.0008 




C 20 


0.013 


0.0085 ± 0.0001 




r + 
c n 


-0.10 


-0.066 ±0.010 




c 01 


0.71 


0.62 ±0.06 


0.51 ±0.10 


c 02 


0.23 


0.34 ±0.03 





Table 3: The subthreshold parameters. The units 
2 ' and m^~ 2j+1 for c£ and c~- respectively. 



are m 



Shown are the results for fit 10 of [45] and the disper- 
sive results from [40] and [28]. 



tribution and shows a reason- 
able convergence over the entire 
region of variation of (L4, Lg) we 
have covered. It is also in good 
agreement with the dispersive cal- 
culation inside the whole region. The result for c^q shows somewhat less good con- 
vergence but it is still acceptable. This is what was used in earlier analyses of irK 
scattering to get a determination of the 1/N C suppressed constant L\. As can be seen 
in Fig. 3(a) the matching with the dispersive values is obtained within the region 
L\ ~ —0.1 10~ 3 with only a fairly weak dependence on the value of Lg. 

The results for the other subthreshold parameters are more difficult to interpret. 
They show a variety of behaviours: i) some subthreshold parameters display indication 
of reasonable convergence while others obviously do not converge, ii) Some agree well 
with the dispersive values while others do so only for large values of L r 4 , Lg. There is 
also no obvious pattern to which values for L£, Lg gives the best agreements. Some of 
these difficulties could be due to the fact that the lowest order is small in the region 
relevant for the subthreshold expansion as can be seen in Figs. 1 and 2. 

We now discuss some of them to show these issues also considering their agree- 
ment with the dispersive calculation. The latter can also be judged from Tab. 3. 
The Cq component agrees with the dispersive result in a very small region for large 
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(a) (b) 

Figure 3: (a) The subthreshold parameter cf as a function of the input values of L\ 
and Lq used in the L\ determination. The unlabeled planes are the result of [40] with their 
errors, (b) The same for the subthreshold parameter Cq . 



negative L\. That was precisely the place where the \ 2 °f the fits for the input 
parameters started getting large, cjj agrees excellently at order p A but gets fairly 
large p 6 corrections. It starts agreeing once more for larger positive values of L\ and 
Lg than considered here. cj~ had a large negative estimate of the contribution from 
the p 6 constants. This drives the total p 6 contribution to be negative and the total 
result stays between 0.04 and 0.15, significantly below the dispersive result of [40]. 
The remaining subthreshold parameters all have large p 6 corrections and it is not 
clear whether we have a convergent series or not. But the general size and the sign 
is correct. 

We now turn to the scattering lengths. The kinematical quantities here have 
values which are already large for ChPT but looking at Figs. 1 and 2 the convergence 
seems fine in that region. The finer features like higher partial waves and the radii 
might however work less well. 

In Fig. 4(a) we have plotted the S wave scattering length in the isospin 1/2 
channel. The series shows a nice convergence and agrees with the dispersive result 
for most of the L\-L r % region we considered. Only a small region of negative L r G 
and positive L\ disagrees. The result for the S wave scattering length in the / = 
3/2 channel, shown in Fig, 4(b), has qualitatively the same behaviour, ruling out a 
somewhat larger region of the (L\, L r G ) plane. For the P wave scattering lengths, we 
get agreement in the 1 = 1/2 channel with the dispersive result in essentially the 
whole region considered with a preference for positive values of L\. The 1 = 3/2 
channel, d^ 2 has large p 6 corrections always leading to a value significantly above 
the dispersive result. Looking at higher threshold parameters the picture is again 
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(a) (b) 

Figure 4: (a) The scattering length Oq 2 as a function of the input values of L\ and L r a 
used in the L\ determination. The unlabeled planes are the result of [40] with their errors, 
(b) The same for the scattering length a . 



mixed, b\j 2 is typically 40 to 60% above the dispersive result, 6q /z is too small by 
20 to 50% and b 1 ^ 2 and b^ 2 have obvious convergence problems. We have shown the 
results for fit 10 of [45] in Tab. 4 together with dispersive estimates of [40]. 

The value of the amplitude at the Cheng- 
Dashen point is related to the kaon sigma 
term. The dependence on L\ and L r 6 is 
shown in Fig. 5 together with the dispersive 
result of [40] . The corrections are large over 
most of the region and can be compared 
with the direct calculation of the sigma term 
shown in Fig. 9(b) of [19]. 

The overall picture we thus obtain in 
the end is rather mixed. If we look only 
at the quantities, cf , Cq , a l J 2 and d^ 2 we 
see a series that converges reasonably well 
and reasonable agreement with the disper- 
sive result from [40] is found for large re- 

Lr) plane. In 



3/2 



Fit 10 



[40] 



1/2 
% 


0.220 


0.224 ±0.022 


in V 2 

10a/ 


0.18 


0.19 ±0.01 


10ao /2 


-0.47 


-0.448 ± 0.077 


10 2 af 


0.31 


0.065 ± 0.044 


106 /2 


1.3 


0.85 ±0.04 


10&o /2 


-0.27 


-0.37 ±0.03 


T+ 

CD 


2.11 


3.90 ± 1.50 



Table 4: The results for the scatter- 
ing lengths and ranges and the amplitude 
at the Cheng-Dashen point as well as the 
dispersive result. The scattering lengths 
and ranges are given in units of m ll +. 



gions of the values in the (L 4 , ±j 6 
particular the value with (L^Lg) ~ (0,0), corresponding to fit 10 of [53], is located 
inside the allowed region. The other quantities present more difficult to interpret 
results depending on how one judges their convergence and the (dis) agreement with 
the dispersive results of [40]. 
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8.2 ttK, 7T7T and Scalar Form-factors 



Even though it seems that the ttK pro- 
cesses alone do not provide strong re- 
strictions on the low-energy constants 
we can use it together with ttti scat- 
tering and the scalar form-factor to re- 
strict the region in the {L\, Lg) plane 
allowed. A full analysis is planned for 
the future but here we discuss the present 
results together with the earlier ones of 
[19, 18]. 

The constraints in the (L r 4 , Lg) plane 
come from several sources: 

1 . The region inferred by the scalar 
form-factor analysis of [19]. 



rK,2r4 




Figure 5: 



K- 0.35 10" 



The amplitude at the Cheng- 
Dashen point as a function of the input val- 
ues of L\ and L\ together with the dispersive 
estimates of [40]. 

This came from two arguments: (i) The assumption that the pion and kaon isoscalar 
scalar form-factors at zero do not deviate by large factors from their lowest order 
values, as judged from Figs. 9(a) and 9(b) in [19]. (ii) The agreement of the ChPT 
calculation of the pion scalar radius with the dispersive results. Notice that the 
dispersive results used the values of the form-factors at zero as input. The ChPT 
prediction for the radius alone was rather constant as shown in Fig. 11(a) in [19]. 

2. From tttt scattering [18] we got constraints from four sources. The parameters 
C\ and C*2 of [12] similar to subthreshold parameters and the scattering lengths Oq 
and ajj. The constraints are (a) a(j agrees reasonably well over the whole region 
considered, see Fig. 4(a) in [18]. (b) ag gives the strongest constraints, see Fig. 4(b) 
in [18]. It requires 



L r 4 > 0.27 10" 



+ 0.9L' 



(8.2) 



(c) C\ imposes L\ > 0.2 10~ 3 and does not constrain Lg, see Fig. 5(a) in [18]. (d) 
C2 does not provide significant more constraints than the above, note that the plot 
shown in [18] is erroneous. 

3. We review now the constraints we found in this article from ttK scattering, 
(a) Cq imposes no constraint, see Fig. 3(b). (b) c^ , delimits 



L\ < 0.28 10" 3 + 0.25 L r 6 , 
1/2 



i.3) 



see Fig. 3(a). (c) The constraint from a^ /z is contained in the a^ 2 , see Fig. 4(a). 
(d) The constraint from aj 2 , see Fig. 4(b), is Lg > L\ — 0.4 10 -3 . 

Here we have used an error double the ones quoted in the dispersive results of 
[12, 40] to take the convergence of the chiral series somewhat into account as well. 
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The nK constraint is essentially only from cf , while for ttty scattering it is 
mainly with a little extra from C\. These constraints do not overlap but the 
various regions almost touch for L r 4 ss 0.2 10~ 3 and L r % < 0.0 10~ 3 . At present the 
situation is only marginally compatible. 

With the uncertainties associated in the calculations: the estimates of 0(p 6 ) 
LECs, correlations between the L\ and the fact that the errors on the L\ from K ei: 
the masses and decay constants is not yet taken into account the above conclusion 
is preliminary. 



9. Conclusions 

In this paper we have calculated the 7rK scattering amplitude to next-to-next-to- 
leading order. We have presented analytically the results to next-to-leading order 
and the dependence on the p 6 LECs C\. The remaining analytical expressions at 
order p 6 are very long and can be obtained from [50] or from the authors. This 
calculation is the main result of this work. 

We presented some numerics with the LECs at order p 4 and p 6 set equal to 
zero at the scale of the rho mass. These results allow a first impression about 
the convergence of the series for various quantities, in particular we have presented 
results for the subthreshold parameters, the scattering lengths and the amplitude at 
the Cheng-Dashen point. 

The second part of this work was a first attempt of extending the order p 4 work 
of constraining low-energy constants from tcK scattering of [40, 39, 58]. In these 
works values for L r 4 were suggested that are positive and different from zero. We 
used in this work as inputs the correlated values for the L[ determined from the 
Ke4 form-factors, pseudoscalar meson masses and decay constants and estimated 
the contributions from the order p e constants with the same procedure as was used 
in previous next-to-next-to-leading processes. In contrast to the p 4 results, a first 
estimate for the quantities which appear to be most reliably obtainable from the nK 
scattering amplitude at order p 6 are fully compatible with both the 1/N C suppressed 
LECs L\ and Lg being equal to zero. Contrary to expectations, the study of these 
nK (sub)threshold parameters alone do not allow to draw definite conclusions on 
the presence of large Zweig rule violating contributions as discussed in [15, 16] and 
references therein. 

The remaining nK quantities present a rather mixed picture, the convergence 
of the series is often questionable and the agreement with the results from the dis- 
persive analysis is at the same level but no obvious large discrepancies exist. The 
estimated contribution of the p 6 constants to many quantities is fairly large and 
rather uncertain, especially those involving the scalars. 
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We have also studied how these results fit together with the earlier ones on txtx 
scattering and the scalar form-factors. We found only marginal compatibility as 
described in Sect. 8.2. 

Planned work for the future is to combine all existing order p 6 calculations in 
three flavour ChPT in order to determine from experiment and/or dispersion theory 
as many as possible of the p 6 constants and to fully take into account all correlations 
for the L\ and errors on the experimental and dispersive inputs used. 

Acknowledgments 

The program FORM 3.0 has been used extensively in these calculations [59]. J.B. 
and P.D. are supported in part by the Swedish Research Council and European 
Union TMR network, Contract No. HPRN-CT-2002-00311 (EURIDICE). PT ac- 
knowledges support by the Spanish Research Council. 

References 

[1] S. Weinberg, Physica A 96 (1979) 327. 
[2] J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142. 
[3] J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465. 
[4] J. Bijnens, Nucl. Phys. B 337 (1990) 635. 

[5] C. Riggenbach, J. Gasser, J. F. Donoghue and B. R. Holstein, Phys. Rev. D 43 
(1991) 127. 

[6] J. Stern, H. Sazdjian and N. H. Fuchs, Phys. Rev. D 47 (1993) 3814 
[arXiv:hep-ph/9301244]. 

[7] J. Bijnens, G. Colangelo, G. Ecker, J. Gasser and M. E. Sainio, Phys. Lett. B 374 

(1996) 210 [arXiv:hep-ph/9511397]. 

[8] J. Bijnens, G. Colangelo, G. Ecker, J. Gasser and M. E. Sainio, Nucl. Phys. B 508 

(1997) 263 [Erratum-ibid. B 517 (1998) 639] [arXiv:hep-ph/9707291]. 

[9] J. Bijnens, G. Colangelo and P. Talavera, JHEP 9805 (1998) 014 
[arXiv:hep-ph/9805389]. 

[10] B. Ananthanarayan, G. Colangelo, J. Gasser and H. Leutwyler, Phys. Rept. 353 
(2001) 207 [arXiv:hep-ph/0005297]. 

[11] G. Colangelo, J. Gasser and H. Leutwyler, Phys. Lett. B 488 (2000) 261 
[arXiv:hep-ph/0007112]. 



- 19 - 



[12] G. Colangelo, J. Gasser and H. Leutwyler, Nucl. Phys. B 603 (2001) 125 
[arXiv:hep-ph/0103088]. 

[13] S. Pislak et al. [BNL-E865 Collaboration], Phys. Rev. Lett. 87 (2001) 221801 
[arXiv:hep-ex/0106071]. 

[14] S. Pislak et al, Phys. Rev. D 67 (2003) 072004 [arXiv:hep-ex/0301040]. 

[15] S. Descotes-Genon, L. Girlanda and J. Stern, Eur. Phys. J. C 27 (2003) 115 
[arXiv:hep-ph/0207337]. 

[16] S. Descotes-Genon, N. H. Fuchs, L. Girlanda and J. Stern, arXiv:hep-ph/0311120. 

[17] L. Girlanda, J. Stern and P. Talavera, Phys. Rev. Lett. 86 (2001) 5858 
[arXiv:hep-ph/0103221]. 

[18] J. Bijnens, P. Dhonte and P. Talavera, JHEP 0401 (2004) 050 
[arXiv:hep-ph/0401039]. 

[19] J. Bijnens and P. Dhonte, JHEP 0310 (2003) 061 [arXiv:hep-ph/0307044]. 

[20] G. Amoros, J. Bijnens and P. Talavera, Nucl. Phys. B 568 (2000) 319 
[arXiv:hep-ph/9907264]. 

[21] P. Post and K. Schilcher, Nucl. Phys. B 599 (2001) 30 [arXiv:hep-ph/0007095]. 

[22] P. Post and K. Schilcher, Eur. Phys. J. C 25 (2002) 427 [arXiv:hep-ph/0112352]. 

[23] J. Bijnens and P. Talavera, JHEP 0203 (2002) 046 [arXiv:hep-ph/0203049]. 

[24] J. Bijnens and P. Talavera, Nucl. Phys. B 669 (2003) 341 [arXiv:hep-ph/0303103]. 

[25] S. Weinberg, Phys. Rev. Lett. 17 (1966) 616. 

[26] R.W. Griffith, Phys. Rev. 176 (1968) 1705. 

[27] C. B. Lang, Fortsch. Phys. 26 (1978) 509. 

[28] C. B. Lang and W. Porod, Phys. Rev. D 21 (1980) 1295. 

[29] N. Johannesson and G. Nilsson, Nuovo Cim. A 43 (1978) 376. 

[30] V. Bernard, N. Kaiser and U. G. Meissner, Nucl. Phys. B 357 (1991) 129. 

[31] V. Bernard, N. Kaiser and U. G. Meissner, Phys. Rev. D 43 (1991) 2757. 

[32] J. Sa Borges, J. Soares Barbosa and V. Oguri, Phys. Lett. B 412 (1997) 389. 

[33] J. Sa Borges and F. R. A. Simao, Phys. Rev. D 53 (1996) 4806. 

[34] V. Bernard, N. Kaiser and U. G. Meissner, Nucl. Phys. B 364 (1991) 283. 



- 20 - 



[35] M. Jamin, J. A. Oiler and A. Pich, Nucl. Phys. B 587 (2000) 331 
[arXiv:hep-ph/0006045]. 

[36] U. G. Meissner and J. A. Oiler, Nucl. Phys. A 679 (2001) 671 
[arXiv:hep-ph/0005253]. 

[37] A. Roessl, Nucl. Phys. B 555 (1999) 507 [arXiv:hep-ph/9904230] . 

[38] M. Frink, B. Kubis and U. G. Meissner, Eur. Phys. J. C 25 (2002) 259 
[arXiv:hep-ph/0203193]. 

[39] B. Ananthanarayan and P. Buttiker, Eur. Phys. J. C 19 (2001) 517 
[arXiv:hep-ph/0012023]. 

[40] P. Buttiker, S. Descotes-Genon and B. Moussallam, arXiv:hep-ph/0310283. 

[41] A. Nehme, Eur. Phys. J. C 23 (2002) 707 [arXiv:hep-ph/0111212]. 

[42] A. Nehme and P. Talavera, Phys. Rev. D 65 (2002) 054023 [arXiv:hep-ph/0107299]. 

[43] B. Kubis and U. G. Meissner, Phys. Lett. B 529 (2002) 69 [arXiv:hep-ph/0112154]. 

[44] B. Kubis and U. G. Meissner, Nucl. Phys. A 699 (2002) 709 [arXiv:hep-ph/0107199]. 

[45] G. Amoros, J. Bijnens and P. Talavera, Nucl. Phys. B 602 (2001) 87 
[arXiv:hep-ph/0101127]. 

[46] A. Pich, Lectures at Les Houches Summer School in Theoretical Physics, Session 68: 
Probing the Standard Model of Particle Interactions, Les Houches, France, 28 Jul - 5 
Sep 1997, [hep-ph/9806303]; 

G. Ecker, Lectures given at Advanced School on Quantum Chromodynamics (QCD 
2000), Benasque, Huesca, Spain, 3-6 Jul 2000, [hep-ph/0011026]; 
S. Scherer, hep-ph/0210398. 

[47] J. Bijnens, G. Colangelo and G. Ecker, JHEP 9902 (1999) 020 
[arXiv:hep-ph/9902437]. 

[48] J. Bijnens, G. Colangelo and G. Ecker, Annals Phys. 280 (2000) 100 
[arXiv:hep-ph/9907333]. 

[49] J. Bijnens, G. Colangelo and G. Ecker, Phys. Lett. B 441 (1998) 437 
[arXiv:hep-ph/9808421]. 

[50] These can be downloaded from http://www.thep.lu.se/~bijnens/chpt.html. 

[51] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321 (1989) 311. 

[52] G. Ecker, J. Gasser, H. Leutwyler, A. Pich and E. de Rafael, Phys. Lett. B 223 
(1989) 425. 



- 21 - 



[53] G. Amoros, J. Bijnens and P. Talavera, Nucl. Phys. B 585 (2000) 293 [Erratum-ibid. 
B 598 (2001) 665] [arXiv:hep-ph/0003258]. 

[54] V. Cirigliano, G. Ecker, H. Neufeld and A. Pich, JHEP 0306 (2003) 012 
[arXiv:hep-ph/0305311]. 

[55] J. Bijnens, G. Colangelo and J. Gasser, Nucl. Phys. B 427 (1994) 427 
[arXiv:hep-ph/9403390]. 

[56] J. Bijnens, E. Gamiz, E. Lipartia and J. Prades, JHEP 0304 (2003) 055 
[arXiv:hep-ph/0304222]. 

[57] M. Knecht and A. Nyffeler, Eur. Phys. J. C 21 (2001) 659 [arXiv:hep-ph/0106034]; 
S. Peris, M. Perrottet and E. de Rafael, JHEP 9805 (1998) 011 
[arXiv:hep-ph/9805442]. 

[58] B. Ananthanarayan, P. Buettiker and B. Moussallam, Eur. Phys. J. C 22 (2001) 133 
[arXiv:hep-ph/0106230]. 

[59] J. A. Vermaseren, math-ph/0010025. 



- 22 - 



